cond-mat/9701159, submitted to Phys. Rev. B 



Theory of a spherical quantum rotors model: 
low— temperature regime and finite— size scaling 

Hassan Chamati 1 , Ekaterina S. Pisanova 2 and Nicholay S. Tonchev 1 ! 

1 Institute of Solid State Physics, Bulgarian Academy of Sciences, 

Tzarigradsko chaussee 72, 1784 Sofia, Bulgaria 

2 University of Plovdiv, Faculty of Physics, 24 "Tzar Assen" str., 

4000 Plovdiv, Bulgaria 
(November 26, 1996) 



Abstract 



The quantum rotors model can be regarded as an effective model for the low- 
temperature behavior of the quantum Heisenberg antiferromagnets. Here, we 
consider a d-dimensional model in the spherical approximation confined to 
a general geometry of the form L d ~ d> xoo d ' xL^. ( L-linear space size and 
L T -temporal size) and subjected to periodic boundary conditions. Due to 
the remarkable opportunity it offers for rigorous study of finite-size effects at 
arbitrary dimensionality this model may play the same role in quantum critical 
phenomena as the popular Berlin-Kac spherical model in classical critical 
phenomena. Close to the zero-temperature quantum critical point, the ideas 
of finite-size scaling are utilized to the fullest extent for studying the critical 
behavior of the model. For different dimensions 1 < d < 3 and < d! < d a 
detailed analysis, in terms of the special functions of classical mathematics, 
for the susceptibility and the equation of state is given. Particular attention 
is paid to the two-dimensional case. 

PACS numbers: 75.10.Jm, 63.70.+h, 05.30-d, 05.20.-y, 02.30 
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I. INTRODUCTION 



In recent years there has been a renewed interest0~i in the theory of zero-temperature 
quantum phase transitions initiated in 1976 by Hertz's quantum dynamic renormalization 
group! for itinerant ferromagnets. Distinctively from temperature driven critical phenomena, 
these phase transitions occur at zero temperature as a function of some non-thermal control 
parameter (or a competition between different parameters describing the basic interaction 
of the system), and the relevant fluctuations are of quantum rather than thermal nature. 

It is well known from the theory of critical phenomena that for the temperature driven 
phase transitions quantum effects are unimportant near critical points with T c > 0. It 
could be expected, however, that at rather small (as compared to characteristic excitation 
in the system) temperature, the leading T dependence of all observables is specified by the 
properties of the zero-temperature critical points, which take place in quantum systems. The 
dimensional crossover rule asserts that the critical singularities of such a quantum system 
at T = with dimensionality d are formally equivalent to those of a classical system with 
dimensionality d+ z (z is the dynamical critical exponent) and critical temperature T c > 0. 
This make it possible to investigate low-temperature effects (considering an effective system 
with d infinite space and z finite time dimensions) in the framework of the theory of finite- 
size scaling (FSS). The idea of this theory has been applied to explore the low-temperature 
regime in quantum systems (see Refs. |5|-[7|), when the properties of the thermodynamic 
observables in the finite-temperature quantum critical region have been the main focus of 
interest. The very quantum critical region was introduced and studied first by Chakravarty 
et ali using the renormalization group methods. The most famous model for discussing 
these properties is the quantum nonlinear 0(n) sigma model (QNLoM). 

Recently an equivalence between the QNLoM in the limit n — > oo and a quantum version 
of the spherical model or more precisely the "spherical quantum rotors" model (SQRM) was 
announced.0 The SQRM is an interesting model in its own. Due to the remarkable oppor- 
tunity it offers for rigorous study of finite-size effects at arbitrary dimensionality SQRM 
may play the same role in quantum critical phenomena as the popular Berlin-Kac spherical 
model in classical critical phenomena. The last one became a touchstone for various scal- 
ing hypotheses and source of new ideas in the general theory of finite-size scaling (see for 
example Refs. |!2}-|T9l and references therein). Let us note that an increasing interest related 
with the spherical approximation (or large n-limit) generatingtractable models in quantum 
critical phenomena has been observed in the last few years.EJ'cHcHl 

In Ref. [11], the critical exponents for the zero-temperature quantum fixed point and the 
finite-temperature classical one as a function of dimensionality was obtained. What remains 
beyond the scope of Ref. |TT] is to study in an exact manner the scaling properties of the 
model in different regions of the phase diagram including the quantum critical region as a 
function of the dimensionality of the system. In the context of the finite-size scaling theory 
both cases: (i) The infinite (^-dimensional quantum system at low-temperatures oo d xL z T 
( 

[L T ~ I J is the finite-size in the imaginary time direction) and (ii) the finite system 

confined to the geometry L d ~ d 'xoo d 'xL z T (L-linear space size) are of crucial interest. 

Earlier a class of exactly solvable lattice models intended to study the displacive 
structural phase transition have been intensively considered in both finite-size and bulk 
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geometry.E2rE3 The main feature of these models is that the real anharmonic interaction 
is substituted with its quantum mean spherical approximation reducing the problem to an 
exactly solvable one. We expect that the analytical technique proposed below will apply to 
these models too. 

In this paper a detailed theory of the scaling properties of the SQRM with nearest- 
neighbor interaction is presented. The plan of the paper is as follows: we start with a brief 
review of the model and the basic equation for the quantum spherical field in the case of 
periodic boundary conditions (Section ||). Since we would like to exploit the ideas of the 
FSS theory, the bulk system in the low-temperature region is considered like an effective 
(d + 1) -dimensional classical system with one finite (temporal) dimension. This is done to 
enable contact to be made with other results based on the spherical type approximation e.g. 
in the framework of the spherical model and the QNLcrM in the limit n — > oo (Section [TIT] ). 
In Section |Iy| we consider FSS form of the spherical field equation for the system confined 
to the general geometry L d ~ d xoc d xL*. This equation turns out to allow for analytic studies 
of the finite-size and low-temperature asymptotes for different d and d! . Special attention 
is laid on the two dimensional system. The remainder of the paper contains the details of 
the calculations: Appendices [A]. [B|. 



II. THE MODEL 



The model we will consider here describes a magnetic ordering due to the interaction of 
quantum spins. This has the following form0 



« = + f E^ 2 - H T, s i 



(2.1) 



where S t are spin operators at site £, the operators V e are "conjugated" momenta (i.e. 
[S e , Sp] = 0, [P e ,V e ,] = 0, and [T-^v&z] = i8 u ,, with % = 1), the coupling constants J e>e , = J 
is between nearest neighbors only,E3 the coupling constant g is introduced so as to measure 
the strength of the quantum fluctuations (below it will be called quantum parameter), H is 
an ordering magnetic field, and finally the spherical field \i is introduced so as to insure the 
constraint 



N. 



(2.2) 



Here N is the total number of the quantum spins located at sites of a hypercubical 
lattice of size L\ x L 2 x ■ ■ • x = N and (■ • •) denotes the standard thermodynamic average 
taken with 7i. 

Let us note that the commutation relations for the operators S £ and V e together with the 
quadratic kinetic term in the Hamiltonian ( |2.1| ) do not describe quantum Heisenberg-Dirac 



spins but quantum rotors as it was pointed out in Ref. [1 1 
Under periodic boundary conditions, Eq. 



takes the form 



A 

2iV 



E 



i 



+ 2Eti(l-cosg i ) 



: coth 



2t\ 



+ 2^(1 - cosgi 



(2.3) 



i=i 
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where we have introduced the notations: A = is the normalized quantum parameter, 
t = 7-the normalized temperature, h = -M= - the normalized magnetic field, b = and 
= j — 2d is the shifted spherical field. 

In Eq. (|2.3| ) the vector q is a collective symbol, which has for Lj odd integers the com- 
ponents: 

f 2im\ 27rn<j 1 f L-l L,- — 1 "l 

{— •-•TTI- nie {-—' ■ '—)■ 

A previous direct analysis0 of Eq. (|2.3| ) in the thermodynamic limit shows that there can 
be no long-range order at finite temperature, for d < 2 (in accordance with the Mermin- 
Wagner theorem). For d > 2 one can find long-range order at finite temperature up to a 
critical temperature t c (X). Here we shall consider the low-temperature region for 1 < d < 3. 



III. THE INFINITE SYSTEM 

In the thermodynamic limit the <i-dimensional sum over the momentum vector q in 
Eq. ( |2.3| ) changes in d integrals over the q^s in the first Brillouin zone and the equation for 
the shifted spherical field <fi reads 



t ~ /■*" r* 1 /i 2 



2 ' 



(3.1) 



After some algebra (see Appendix [S]), Eq. ( |3.1| ) takes the form (1 < d < 3) 



A A c (47r)( rf + 1 )/ 2 



r 



1 - rf N 



,(d-l)/2 



(47T 



.(d+l)/2 



2 ' 2t 



(3.2) 



where A c is the quantum critical point and 

' d- 1 



l,y) =2Y,(y™y u K v (2my). 



m=l 



(3.3) 



Here K u (x) is the MacDonald function (second modified Bessel function). The asymptotic 
forms of the functions /Ci 1, are studied in Appendix [FJ. In the remainder of this 
section we will study the effect of the temperature on the susceptibility and the equation of 
state near the quantum critical fixed point. 



A. Zero field susceptibility 

After making vanishes the field h, from Eq. ( p.2|) we find that the normalized zero-field 
susceptibility x — _1 on the line A = A c (t — > + ) is given by 
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X = ht~'\ (3.4) 



4y, 
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where yo is the universal solution of 



d 



2K 



'd-l 



(3.5) 



The behavior of the universal constant yo as a function of the dimensionality d of the system 
is shown in FIG. [1]. 

One can see that the low-temperature behavior of the susceptibility increases as the 
inverse of the square of the temperature above the quantum critical point. 

In what follows we will try to investigate Eq. (|3.2| ) for different dimensions (1 < d < 3) 
of the system and in different regions of the (t, A) phase diagram. Introducing the "shifted" 
critical value of the quantum parameter due to the temperature by 



1 



i 

A^ T 2vr( d + 1 )/ 2 VA C 



d-i 



1 r 



'd-l' 



\C(d-l)\, 



(3.6) 



(where ((x) is the Rieman zeta function) one has to make a difference between the two cases 
d < 2 "sign — " and d > 2 "sign +". In the first case (1 < d < 2), it is possible to define the 
quantum critical region by the inequality 



1 



1 



d-l 



27r (d+l)/2 \\ c (t) 



r 



i 



|C(d-l)|. 



(3.7) 



For 1 < d < 2 the function K{y,y) ~ y 1 and by substitution in (|3.2|) we obtain for 
A < A c (outside of the quantum critical region) 



X 



r(i 



(A7T) d / 2 SX 



2/(d-2) 



f 2/(d-2) 



(3.8) 



where 



5\ 



1 1 

a; - a - 

In Eq. ( |3.8|) we see that the susceptibility is going to infinity with power law degree when the 
quantum fluctuations become important (t — > + ) and there is no phase transition driven 
by A in the system for dimensions between 1 and 2. 
In the second case (2 < d < 3), one has 



X 



r(i 



A c (t) 



(4yr) rf / 2 A-A c (t) 



2/(d-2) 



f 2/(d-2) 



(3.9) 



as a solution for A less than A c and greater than the critical value A c (t) of the quantum 
parameter. Here for finite temperatures there is a phase transition driven by the quantum 
parameter A with critical exponent of the <i-dimensional classical spherical model 7 = -Aj. 
This however is valid only for very close values of A to \(t). For A < A c (t) the susceptibility 
is infinite. 

In the region where A > A c the zero-field susceptibility is given by 
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X 



(4tt) 

r 



(d+l)/2 



l-d 



-5X 



2 
l-d 



(3.10) 



This results is valid for every d between the lower and the upper quantum critical dimensions 
i.e. 1 < d < 3. 

The important case d = 2 can be solved easily and one gets 



.V^arcsinhUexp 
A 2 



2ttA 



sx 



For the susceptibility, Eq. ( |3.11D yields 



4vrA, 



X ^ — exp 



SX 



(3.11) 



(3.12a) 



for f 
3.1114. 



A- 1 



^> 1 and A < A c i.e. in the renormalized classical region. For A = A c 



x e 2 \t. 



(3.12b) 



where the universal constant 




.962424.. 



(3.12c) 



was first obtained in the framework o: 
one finite dimension.^ Finally for ^ 
region. 



the 3-dimensional classical mean spherical model with 
3> 1 and A > A c i.e. in the quantum disordered 



X ~ [4:71 5 X] 



-2 



1 + 



2t 



nX c 6X 



exp 



4ttA c 



5X 



(3.12d) 



The first term of Eq. ( |3.12d| ) is a particular case of Eq. Q3.10| ) for d = 2. 

From Eqs. ( p,12|) one can transparently see the different behaviors of x{T) in tree regions: 
a) renormalized classical region with exponentially divergence as T — >• 0, b) quantum critical 
region with x{T) ~ T~ 2 and crossover lines T ~ |A — A c |, and c) quantum disordered 
region with temperature independent susceptibility (up to exponentially small corrections) 
as T —> 0. The above results ( |3.11| ) and (ft.!2|) coincide in form with those obtained in 
Refs. for the two dimensional QNLaM in the n —>■ oo limit. The only differences are: in 
Eq. (|3.11| ) the temperature is scaled by A, and the critical value A c is given by Eq. ( |A13| ), 
while for the QNLaM it depends upon the regularization scheme. 
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B. Equation of state 

The equation of state of the model hamiltonian Q2.1| ) near the quantum critical point is 
obtained after substituting the shifted spherical field <f> by the magnetization M. through 
the relation 



M 



(3.13) 



in Eq. (|3.2|), which allows us to write the equation of state in a scaling form 



SX 



+ (4tt 



-(d+l)/2 



h 

M 1 



1 1/7 



r 



1 -d 



2/C 



d-1 X ( M v/(3 " 



2 ' 2 



t 



'_h_ 
M 1 



1. 



(3.14) 



We conclude that near the quantum critical point Eq. ( |3.14j ) may be written in general forms 



as 



M s f h (SXM- 1/l3 , (t/xy^M- 1 ^) 



(3.15a) 



or 



M 



X 



f M (5XM- 1/f3 ,hM- 5 ) 



(3.15b) 



In Eqs. ( |3.15| ) fh{x,y) and /x(x,?/) are some scaling functions, furthermore 7 = ^y, 
v = -tzTj P — \ an d 5 = ^1 are ^ ne familiar bulk critical exponents for the (d + 1)- 
dimensional classical spherical model. Eqs. ( |3.15| ) are direct verification of FSS hypothesis 
in conjunction with classical to quantum critical dimensional crossover. They can be easily 



transformed into the scaling form (Eq. (21)) obtained in Ref. |Tl|, however here they are 
verified for 1 < d < 3 instead of 2 < d < 3 (c.f. Ref. |TT1) , i.e. the non-critical case is 
included. 

Hereafter we will try to give an explicit expression of the scaling function fh(x,y) (x = 
5XA4~ 2 , y = (t/A) d_1 A4~ 2 ) in the neighborhood of the quantum critical fixed point. This 



may be performed, in the case j^h/M. -C 1, with the use of the asymptotic form of /C(z/, y) 
the following result for the scaling function 



fh(x,y) 



(4tt 



,d/2 



y 



1 + X + 



1 



2/(d-2) 



r (1 - |J" V 2vr( (i + 1 )/ 2 ~ V 2 

For the special case d = 2 the scaling function is given by the expression 



(3.16) 



h(x,y) = 4y 2 



arcsinh-exp f 2n — — - 



(3.17) 



At x = and ?/ 3> 1 (fixed low-temperature and h — > + ) Eq. ( p,17| ) reduces to 
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/a (0,J/) ~ y 2 exp 



Air 



(3.18) 



In the region x < — 1, and for ?/<l (fixed weak field and t — > + ) the corresponding scaling 
function is 



x + 1 



(3.19) 



and in region x > — 1 and ?/<! we have 



f h {x,y) « 16vr 2 (x + l) 



7r(l + X 



exp — 47r 



1 + x 



(3.20) 



This identifies the zero temperature (y = 0) form of the scaling function ( |3.20D with those 
of the 3-dimensional classical spherical model. 



IV. SYSTEM CONFINED TO A FINITE GEOMETRY 

When the model Hamiltonian ( |2.1| ) is confined to the general geometry L d ~ d 'xoo d 'xL T , 



with < d! < d, equation (2J3) of the spherical field (p takes the form (for a derivational 
details see Appendix |A]) 



1 = 1- (47T)-( d+1 )/ 2 

A \ r 



r 



1 -d 



<p( d -^ > Ky-w [0 1/2 {(Am/t) 2 + (L| Z|) 2 } 1/2 

(27r )(d+l)/2 m ^_ d/) r {(AmA)2 + (L | Z | )2} l/2 



+ 



II 2 



(d-l)/2 1 02' 



(4.1) 



where 



|i| = + £ + • •• + £_* 



1/2 



and the primed summation indicates that the vector with components m = l\ = l 2 
ld-d' = is excluded. 



A. Shift of the critical quantum parameter 



The FSS theory (for a review see Ref. ^8]) asserts, for the temperature driven phase 
transition, that the phase transition occurring in the system at the thermodynamic limit 
persists, if the dimension d' of infinite sizes is greater than the lower critical dimension of 
the system. In this case the value of the critical temperature T c (oo) at which some thermo- 
dynamic functions exhibit a singularity is shifted to T C (L) critical temperature for a system 
confined to the general geometry L d ~ d xoo d , when the system is infinite in d' dimensions and 
finite in [d — d') -dimensions. In the case when the number of infinite dimensions is less than 
the lower critical dimension, there is no phase transition in the system and the singularities 
of the thermodynamic functions are altered. The critical temperature T c (oo) in this case is 
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shifted to a pseudocritical temperature, corresponding to the center of the rounding of the 
singularities of the thermodynamic functions, holding in the thermodynamic limit. 

In our quantum case, having in mind that we have considered the low-temperature 
behavior of model (|2.1|) in the context of the FSS theory it is convenient to choose the 
quantum parameter A like a critical instead of the temperature t and to consider our system 
confined to the geometry L d ~ d xoo d xL T . So the shifted critical quantum parameter X c (t, L) = 
XtL is obtained by setting <fr = in Eq. This gives 



j_ i r(¥) 

XtL X c 47r( d +i)/2 



E ((X tL m/tr + (L\l\ff- d)/2 . (4.2) 

m,l(d-d') 



The sum in the r.h.s. of Eq. ( |4.2| ) is convergent for d' > 2, however it can be expressed 
in terms of the Epstein zeta function 




L 2 Z 2 +^Vm 2 ;d-l| = £ \L 2 l 2 +\'-) m 2 \ , (4.3) 



which can be regarded as the generalized (d — d' + l)-dimensional analog of the Riemann 
zeta function C(^rO ( see Ref- P§D - In the case under consideration the Epstein zeta function 



has only a simple pole at d' = 2 and may be analytically continued for < d' < 2 to give 
a meaning to Eq. ( |4.2| ) for d' < 2 as well. It is hard to investigate the sum appearing in 

2 72 I I r2/2 i /A\ 2 „2 



Eq. ( |4.3| ) . The anisotropy of the sum L 2 lf + ■ • ■ + L 2 l 2 d _ d , + lj j m 2 is an additional problem. 
That is for what we will try to solve it asymptotically, considering different regimes of 
the temperature, depending on that wether L <C ^ or L > ^ which will be called, 
respectively, very low-temperature regime and low-temperature regime. 



1. Low-temperature regime <C L 
In this case after some algebra the resulting expression is 

, f -sd/2 rl-d/2 / oo / '\ d / 2 - 1 / + x 

+ f ) E E 2T f Lm|i • (44) 

The first term of the rhs of Eq. (|4.4| ) is the shift of the critical quantum parameter (see 
Eq. ( P-6| )) due to the presence of the quantum effects in the system. The second term is 
a correction resulting from the finite sizes. It is just the shift due the finite-size effects in 
the o?-dimensional spherical mo multiplied by the temperature scaled to the quantum 
parameter. Here the (d — d')— fold sum may be continued analytically beyond its domain 
of convergence with respect to "d" and "cf" (which is 2 < d' < d). The last term is 
exponentially small in the considered limit i.e. ^ <^ L. 

In the borderline case d = 2, Eq. ( |4.4| ) reduces to 
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A 



tL 



\ r 2vrA 



t I . *L 

5 + 7i? + In 



fL 



2A 



2 E E^o 2tt— Lm\l 

IL l(2-d')m=l ^ *tL 



t 



(4.5) 



where ^ E = .577... is the Euler constant and B is a constant depending on the dimensionality 
d-d'. 

In the particular case of strip geometry (d! = 1), B = — ln27r and the shift is given by 



1 



1 



t 



A c 2ttA 



tL 



[lE + In 



tL 



4vrA tL 



For the fully finite geometry case (d! = 0), B 



In 



A 



tL 



X r 2vrA 



t I , 

7e + In ■ 



4vrA 



In 



2v^ 



(4.6) 



and we obtain for the shift 



tL 



47Ty / 7r 



(4.7) 



Let us note that in the rhs of Eqs. ( |4.6| ) and ( f4.7|) exponentially small corrections are omitted. 



2. Very low-temperature regime 4 3> £ 



The final result in this case is given by the expression 



1 

A*L 



1 

AT 



47r (d+l)/2 



E l« 

I(d-cZ') 



L d'-d 



27r (d'+l)/2- 



A 



fL 



d'-l 



C(rf' - 1) 



L l/2+d'/2-d , t x ^ 



7T 



d'-l 

T~~ ] E' E (-) 2 (2^m|Z| 



(4.8) 



Here, in the rhs, the first term is the expression of the shift of the critical quantum pa- 
rameter, at zero temperature ,0 due to the finite sizes of the system. This is equivalent to 
the shift of a (d + l)-dimensional spherical model confined to the geometry L d + l ~ d xoo d . 
The second term gives a corrections due to the quantum effects. This is the shift of crit- 
ical quantum parameter of a (i'-dimensional infinite system multiplied by the volume of a 
(d — d') -dimensional hypercube. The third term is exponentially small in the limit of very 
low temperatures. For d' — 1 Eq. (|4.8|) yields 



1 

AtL 



2n 



where 



Co 



27r (<*-l)/2 



r 



l(d-l) m=l 

(d-l 



(4.9) 



and the expressions for the constants C are quite complicated expect for some special cases 
(see Refs. e.g. for d = 2, d! = 1 -e we have C = 7 B -ln(47r) (c.f. Ref. Eq. (30.104)). 
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In the particular case d! = 1 and d = 2 the constant Co is given by Co = 7^ — In Arc and 



a!z-^SiK +1d ^}' (410) 



Comparing between Eqs. (O) and ( 4.10 ) one can see the crucial role (in symmetric form) 



of L or \th/t in the low-temperature regime and the very low-temperature one, respectively 
In the other particular case of a two-dimensional bloc geometry d! = and d = 2, from 
Eq. Q4.8| ) one gets (up to exponentially small corrections) 

where 

000 = E 



^(2Z + 1)*- 

Instead of the previous case of low-temperature regime here the lower quantum critical 
dimension d! — 1 is responsible for the logarithmic dependence in Eq. (|4.9|) . This is the 
reason for the significant difference between Eqs. ( |4.10| ) and ( [4.11| ). 

The obtained equations for \ t i will be exploited later for the study of the two dimensional 
case. 



B. Zero field susceptibility 



It is possible to transform Eq. ( |4.1|) in the following equivalent forms 



where 



\ L 2 <p ) ~ (4vr)( d + 1 )/ 2 



1 -d" 



2/C 



'd-l 



tL \ 2 



d-d' + 1 



L(j) 



1/2 < 



(4.12a) 



/ K u (2ysjl 2 + a 2 m 2 

JC a (u\ Pl y)= ]T / ,„ y A V>0, l 2 = ll + ll + --- + lU- (4.12b) 



*(p-i) [yJl 2 + a 2 m 2 



or 



l-d 



5X 



h /t\ 2 fx 



(d+3)/2 



<f)\\J \t 



{Xcj) 1 / 2 /ty- 1 

(4vr)(rf+i)/2 

'd-l 



r 



1 - cf 



— 2/C tL 

A 



d — d' + 1, 



A0 



1/2' 



2f 



(4.13a) 



where 
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fc a (v\p,y) = £i(Hp; a v) 



/ K v [2y\Ja 2 l 2 + m 2 



m,i(p-i) ?/ y a 



(4.13b) 



The functions K a {v\p\y) and /C a (z/|p;j/) are anisotropic generalizations of the /C-function 
introduced in Ref. 0. 

Eqs. ( 4.12 ) and ( [4.13 ) show that the correlation length £ = <\>~ l l 2 will scale like 



or like 



A 
1 



a J 'T : 



A, 



(4.14a) 



(4.14b) 



which suggests also that there will be some kind of interplay (competition) between the 
finite-size and the quantum effects. 

Hereafter we will try to find the behavior of the susceptibility \ = as a function 
of the temperature t and the size L of the system. For simplicity, in the remainder of this 
section, we will investigate the free field case (h — 0). 

1) For \^l 2 < 1, after using the asymptotic form of the function defined in (|4.12b|) (see 
Appendix g) Eq. (Q reads {d! ^ 2, 1 < d < 3) 



+ 



47r (d+l)/2 



r 



E 

m,I(d-d') 



^m) +(LZ) 2 



(4.15) 



Now we will examine Eq. ( |4.15[ ) in different regimes of t and L and for different geometries 
of the lattice: 

a. jcj) 1 / 2 <C 1 and ^ ^> 1: In this case Eq. ( [4.15[ ) transforms into (up to a exponentially 
small corrections c.f. with Eq. ( |4.4| )) 



f jd'-d 

o = 5\ + - . „,„ ru 



A(47T) d '/2- 

A 4^7^ \2 



d'-2)/2 



+ 



1 



27r (rf+l)/2 V A 



(4.16) 



This equation has different type of solutions depending on that whether the dimensionality 
d is above or below the classical critical dimension 2. 

At A = A c and when a" < 2 < d < 3 (i.e. when there is no phase transition in the system) 
we obtain for the zero-field susceptibility 



X 



t \ - 2 /tL\ 2{ - d ~ d ' )/{ - 2 ~ dl) 

\J VaTJ 



r 



d-l 



71 



(d-d'+l)/2 



ni-f 



2 

2-d' 



(4.17) 
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However for 1 < d < 2, Eq. ([4.16 ) has no solution at A = A c obeying the initial condition 
^0 1 / 2 «1. 

Eq. ( 4.17 ) generalizes the bulk result ( |3.4j ) for d close to the upper quantum critical 
dimension i.e. d = 3. 

At the shifted critical quantum parameter A c (t) given by Eq. ( |3.6| ) we get 



X = L 2 



id'-2 



r(|-i) 



7^ £ I' 



2-d 



7r (d-d')/2 F d> 



2 

2-d' 



(4.18) 



However this solution is valid only for 3 > d > 2 > d' i.e. here again there is no phase 
transition in the system. 

b. j(p 1/2 < 1 and ^ < 1: In this case, Eq. ( |4.15|) gives (up to exponentially small 
corrections c.f. with Eq. (|4.8|) ) 



„ t L d '- d „ / d r 
= 6X + - ,. ; „,„ r I 1 - 



A (4vr) d '/2 



2)/2 



l-d 



+ 



d- 1 



4 7r (d+l)/2 



E I' 

J(d-d') 



l-d 



rd'-d 

+— r 



'<£ - l\ ft 



d'-i 



C(d'-i). 



(4.19) 



Here we find that the solutions of Eq. (|4.19|) depend upon that whether the dimensionality 
d' < 1 or d' > 1. 

At A = A c and for 1< d' < 2, Eq. ((D]§ has 



2/(2-d') 



id'-2 



^(d— d'+l)/2 p / 1 ,p\ \ l 



il— d 



2 

2-d' 



(4.20) 



as a solution. For < d' < 1, however, it has no solution obeying the initially imposed 
restriction ^fcf) 1 / 2 1. 

At the shifted critical quantum parameter A C (L) given by§l 



l-d 



A \ C (L) 47r( d + 1 )/ 2 



r 



d- 1 



E \i 

i(d-d') 



l-d 



Eq. ( [4. 19|) has a solution obeying the initial condition ^0 1//2 <C 1 only for d' = 1 + £ and in 
this case the susceptibility behaves like 



X 



1 X] 
(tts) 2 t 2 



1 - e ( 7 E + In - 



(4.21) 



2) For Lcj) 1 ! 2 <C 1, from Eqs. ( |4.12j ) and Eq. ( |B9| ) we get once again Eq. ( [4.16Q . In spite 
of the fact that we have the same equation as in the case j^ 1 ^ 2 <C 1, the expected solutions 
for the susceptibility may be different because of the new imposed condition. Here also we 
will consider the two limiting cases of low-temperature and very low-temperature regimes. 



13 



a. <C 1 and ^ 1: In this case Eq. (|4.15|) again is transformed into Eq. (|4.16|) and 

we obtain at A = A c the solution given by Eq. ( f4.17| ), which is valid only for d' < 2 < d < 3, 
i.e. we have the same solution as in the previous case i.e. jcp 1 ^ 2 <C 1. 

At A = A c (£), we formally obtain Eq. ( 4.18 ) which, however, may be considered like a 
solution only in the neighborhood of the lower classical critical dimension d = 2. For the 
cylindric geometry (d' = 1 and d = 2 + e) we get 



X 



ln47r) 



1 2 



(4.22) 



This result is contained in Eq. (30.109) of Ref. |3l] in the large n-limit case for the NLQcxM. 
In the case of slab geometry d — d' = 1 (d = 2 + e, d' = 1 + s) instead of ( f4.22| ) we obtain 



X 



[1 -e{i E - In 2) -elnef . 



(4.23) 



In the case of a bloc geometry (d = 2 + e and d' = 0) we find the following behavior for 
the susceptibility 



X 



2ne 



( 



1e 



In 



V 



47T 2 



For the case of " quasi-bloc geometry" (d — 2 + e and d' = e) we get 



(4.24) 



X 



2ne 



( 



1 - 



2 lE + In 



2 In 



V 



2^ 



(4.25) 



Here the appearance of e in the denominator in formulas ( 4.21 )-( 4.25 ) signalizes that the 
scaling in its simple form will fail at e — 0. 

b. L(f) 1/2 < 1 and y < 1: Here we find that Eq. (p~T9"l) is valid, and it has Eq. ( ^20]) as 
a solution at A = A c and for < dl < 1. For 1 < cf < 2 the susceptibility is given by 



X \2t. 



ttV2 r (l - f ) 



fvC(rf' - 1) 



2 

2-d' 



(4.26) 



At the shifted critical point A C (L), for the susceptibility we obtain Eq. (|4.26| ) under the 
restriction 2 > d' > 1, which guarantees the positiveness of the quantity under brackets. 

When A < A c for 1 < d < 3 and d' < 2, i.e. when there is no phase transition in the 
system, we obtain 



X 



(47T 



r i 



2 

2-d' 



£-2/ (2-d') j2{d-d')l{2-d') 



(4.27) 
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If d' > 2 there is a phase transition in the system at the shifted value of the critical 
quantum parameter \ t i (the shift in this case is due to the quantum and finite-size effects) 
and Eq. ( [4.1 5| ) transforms to 



_ L d '- d 
1 " A^ ~ t (4^ 



r (i - d - ) <p^\ 



which has the following solutions 

(4tt// 2 



X 



(4n) d / 2 /, X \ 2 - d t -2/{2-d>) L 2(d-d')/(2-d') x > JV 

r(i-f) V x tL )\ 1 h ' 

o, A < A tL 



(4.28) 



(4.29) 



Let us notice that Eqs. fl4.27| ) and ( (4.29| ) are the finite-size forms, for the susceptibility, of 
Eqs. ( |3.8j ) and ( |3.9| ), respectively, found for the bulk system. 



C. Two dimensional case 

The two-dimensional case needs special treatment because of its physical reasonability 
and the increasing interest in the context of the quantum critical phenomena.i~0 From 



Eq. (4T) for d = 2 and in the absence of a magnetic field h = we get 

1/21 



6 i/2 i / exp 

5X = ^ V 



.01/2 ^ m 2 + L 2f 



l7r ^mm-d>) (£ m i + m^ 1/2 



(4.30) 



Introducing the scaling functions Yf = j4> 1 ^ 2 and Yf = L<p 1 ^ 2 , where the superscript d' 
denotes the number of infinite dimensions in the system, and the scaling variable a — ^ it 



is easy to write Eq. ( |4.30| ) in the scaling forms given in Eqs. ( |4.12| ) and (|4.13|) . The solutions 



of the obtained scaling equations will depend on the number of the infinite dimensions in 
the system. Here we will consider the two most important particular cases: strip geometry 
d' = 1 and bloc geometry d' = 0. Our analysis will be confined to the study of the behavior 
of the scaling functions at the critical value of the quantum parameter A c , and at the shifted 
critical quantum parameter X tL (see Section [IV A|) . It is difficult to solve Eq ( |4.30|) by using 



an analytic approach that is for what we will give a numerical treatment of the problem. It 
is, however, possible to consider the two limits: a> 1 i.e. the low-temperature regime and 
a C 1 i.e. the very low-temperature regime. 

Strip geometry (d' = 1): In this case in the rhs of Eq. (|4.30| ) we have a two-fold sum 
which permits a numerical analysis of the geometry under consideration. FIG. |2] graphs the 
variation of the scaling functions and Y£ against the variable a at A = A c . This shows 
that for comparatively small value of the scaling variable a ~ 5 the finite-size behavior 
(see the curve of the function Y^(a)) merges in the low temperature bulk one, while the 
behavior of Y^ia) shows that for relatively not very low-temperatures (a ~ |, L- fixed) the 
system simulates the behavior of a three-dimensional classical spherical model with one 
finite dimensions. The mathematical reasons for this are the exponentially small values of 
the corrections, as we will show bellow. 
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Bloc geometry (d! = 0): In this case the three-fold sum in the rhs of Eq. (|4.30|) is not 
an obstacle to analyze it numerically. For A = A c the behavior of the scaling functions 
Y®(a) and Y t °(a) is presented in FIG. H They have the same qualitative behavior as in the 
strip geometry the only difference is the appearance of a new universal number for t = 0, 
i.e. Q, instead of the constant B as a consequence of the asymmetry of the sum in the 
low-temperature and the very low-temperature regimes. 

Analytically for arbitrary values of the number of infinite dimensions d! we can treat 
first the problem in the low-temperature regime (a ^> 1). In this limit Eq. (|4.30| ) can be 
transformed into (up to small corrections 0(e~ 27Ta )) 

5A = -i-ln2sinhA i/ 2 _ JL £ Kq (l^ 2 \1\) . (4.31) 



For A = A c Eq. (|4.31|) has the solution 

^ i/2 ^ e+ ( 2 - d, )/i(^) 1/2 ^(-^ e ) <« 2 ) 

i.e. the finite-size corrections to the bulk behavior are exponentially small. 

In the very low-temperature regime (a <C 1), Eq. (|4.30|) reads (up to 0(e~ 2w ^ a )) 

<b 1/2 L- 1 ' exp i-L(t) ll2 \l\] L d '- 2 ( A\ ^ 00 (\ , ,„\ 

M = 17 " 17 J Sl + W) £MH (4 - 33) 

which has the solutions: 



for cf = 1, and 

A L L 

for <f = 0. Here Q = 1.511955... is an universal constant. 

In Section [IV A| an analytic continuation of the shift of the critical quantum parameter 



for d = 2 was presented. It is possible to consider the solutions of Eq. ( 4.30|) at A = X t L 
(from Eqs. Q4.6| ), (|4.7|) , ( |4.10|) and (|4.11| )) and for different geometries. In this case the 
scaling functions Y^, Y£, Y t ° and Y® are graphed in FIG. |3|. For df — 1 again we see that a 
symmetry between the two limits aCl and a ^> 1 take place, since the scaling functions 
Y^ and Y£ are limited by the universal constant H. The asymmetric case dl — has two 
different constants S t and limiting the solutions of Yj and from above. 

The constants 2, E t and £l are obtained from the asymptotic analysis (with respect to 
a) of Eq. (|O0| ) for A = A tL . 



In the limit a ^> 1 for arbitrary values of cf we get (from Eq. (|4.31|) ) 
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B + j E + In 



1/2 



E Ko(L<j> 1/2 \l 



(4.36) 



1(2- d') 



where the equation of X t L from Eq. (|4.5|) is used. Eq. ([4.361) has the solutions: 



5 for d! = 1, 
S L for tf = 0, 



(4.37) 



where the universal numbers 



= 7.061132... and = 4.317795... are the solutions of the 
scaling equation ( |4.36| ) for d! — 1 and ef = 0, respectively. 

In the opposite limit a -C 1, for c?' = 1, we get from Eqs. (|4.10 ) and Q4.33Q the equation 



7s + ln 



4vrt 



2 £ Ko (^<P 1/2 m) 

771=1 V 1 / 



which has 



A 



t& -1/2 



as an universal solution. For cf = we have 



A 



^V/ 2 - 6 ) exp f ^V/ 2 ) - _ 6 = o 



A, 



exp 



(4.38) 



(4.39) 



(4.40) 



obtained from Eqs. Q4.11J) and ( |4.33| ), where we have used the identity (P 1 1|) . 
From Eq. (4.40) we obtain the universal result 



r^ v -i/2 
t x 



Et = 6.028966. 



(4.41) 



We finally conclude that if we take A = A c the scaling functions Yf and Yf have similar 
qualitative behavior weakly depending on the geometry (i.e. bloc d! = or strip d! = 1) 
of the system. However, for a given geometry one distinguishes quite different quantitative 
behavior of the scaling functions depending on that whether the quantum parameter A is 
fixed at its critical value, i.e. A = A c , or takes "running" values \l obtained from the "shift 
equations" 0, Q, ( gig ) or (PI|) . 



D. Equation of state 



The equation of state of the model hamiltonian ( |2.1| ) for dimensionalities 1 < d < 3 is 
given by (see Eqs. ( |3.13|) and (|4.1| )) 



= 5\ - (4vr)-( d+1 )/ 2 

h x^- 1 )/ 2 



(27r)( d+1 )/ 2 



E 



' K(d-l)/2 



te) 1/2 ((AmA) 2 + (L|/|) 2 ) 1 



/2 



m,l(d-d') 



1/2 



((Am/t^^l) 2 ) 1 / 2 



(d-l)/2 



(4.42) 
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It is straightforward to write this equation in a similar form as in Eq. Q4.12|) or Eq. fl4.13|) , 
i.e. 



or 



h = M s fli [ 5\M-^, t ±^ 1/u M -W\ 



(4.43a) 



(4.43b) 



Eqs. ( |4.43|) are generalizations of Eqs. (|3.15| ) in the case of systems confined to a finite 
geometry. The appearance of an additional variable y is a consequence of the fact that the 
system under consideration may be regarded as an "hyperparallelepiped" (in not necessary 
an Euclidean space) of linear size L in d — d! directions and of linear size L T in one direction 
with periodic boundary conditions. 



V. SUMMARY 

Since exact solvability is a rare event in statistical physics, the model under consider- 
ation yields a conspicuous possibility to investigate the interplay of quantum and classical 
fluctuations as a function of the dimensionality and the geometry of the system in an exact 
manner. Its relation with the QNLcrM in the n — > oo limit may serve as an illustration of 
Stanley's arguments of the relevance of the spherical approximations in the quantum case. 
Let us note, however, that the use of such type arguments needs an additional more subtle 
treatment in the finite size case. For this reason we gave a brief overview of the bulk low- 
temperature properties, which are similar to those obtained by saddle-point calculation for 
the QNLaM (see Section [II]). 

The Hamiltonian ( |2.1[ ) can be obtained also in the "hard-coupling limit" from a more 
realistic model that take into account quartic self-interaction term Qj, by the anzatz Qg =>■ 
J2e Qj frequently used in the theory of structural phase transitions (see Refs. p3TpB|) . 
So the Hamiltonian ( |2.1|) can be thought as a simple but rather general model to test 
some analytical and numerical techniques in the theory of magnetic and structural phase 
transitions. The discussion of the obtained, in Section |T|, results serves as a basis for 
the further FSS investigations. Identifying the temperature, which governs the crossover 
between the classical and the quantum fluctuations as an additional temporal dimension 
one makes possible the use of the methods of FSS theory in a very effective way. 

In Subsection [IV A| the shift of the critical quantum parameter A as a consequence of the 



quantum and finite-size effects is obtained. In comparison with the classical case (for details 
see Ref. and references therein) here the problem is rather complicated by the presence of 
two finite characteristic lengths L and L T . We observe a competition between finite size and 
quantum effects which reflects the appearance of two regimes: low-temperature (L ^> L T ) 
and very low-temperature (L <C L T ). The behavior of the shift is analyzed in some actual 
cases of concrete geometries e.g. strip and bloc. 

In the parameter space (temperature t and quantum parameter A), where quantum 
zero point fluctuations are relevant, there are three distinct regions named " renormalized 
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classical", "quantum critical", and "quantum disordered". The existence of these regions 
in conjunction with both regimes: low-temperature and very low-temperature make the 
model a useful tool for the exploration of the qualitative behavior of an important class of 
systems. 

In Subsection [IV B| the susceptibility (or the correlation length) is calculated and the 



critical behavior of the system in different regimes and geometries is analyzed. We have 
studied the model fl2.ip via ^-expansion in order to illustrate the effects of the dimensionality 
d on the existence and properties of the ordered phase. An indicative example is given by 
Eqs. (4.22) and ( |4.23 ), while the former is known (see Ref. [H]) the last one is quite different 



and new. These shows that one must be accurate in taking the limit e — » + in the context 
of the FSS theory. 

In Subsection [IV G , special attention is paid to the two-dimensional case. The two 
important cases of strip and bloc geometries are considered. The universal constant G given 
by Eq. ( |3.12c| ), which characterizes the bulk system, is changed to a set of new universal 
constants: (see Eq. S and S L (see Eq. (|Q7D ), and E t (see Eq. flOID). The 



appearance of new universal constants reflects the new situation, when there are two relevant 
values of the quantum parameter A: A = A c in the bulk case and A = Xn in the case of 
finite geometries. Due to their universality these constants may play an important role even 
in studying more complicate model Hamiltonians. The behaviors of the scaling functions 
at the bulk critical quantum parameter A c and the shifted critical quantum parameter Atr, 
are given in FIGs. |2] and |3|. L T is the main characteristic length and the \ corrections are 
exponentially small in the case of low-temperature regime, and vice versa in the case of very 
low-temperature regime. 

The equation of state, for the system confined to the general geometry L d ~ d xoo d xL T , 
is obtained in Subsection |1V O. This reflects the modifications of the scaling functions as a 



consequence of the finite sizes and the temperature. 

Finally let us note that this treatment is not restricted to the model Hamiltonian ( |2.1| ) . 
but it can be applied to a wide class of finite lattice models (e.g. directly to the anharmonic 
crystal model see Refs. ^3|-[26|) and it can also provide a methodology for seeking different 
quantum finite-size effects in such systems. 
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APPENDIX A: 



In this appendix we will derive Eqs. (|3.2|) and (|4.1| ) of the shifted spherical field for 
the model Hamiltonian (|2.1| ) confined to the general geometry L d ~ d xoc d xL T , with periodic 
boundary conditions, in the low temperature regime. To achieve that, let us start with 
Eq. (O) 
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l = W d (0,L,t) + ^, 



(Ala) 



where we have used the notation 

t 00 



w d {<j>,L,t) = - y: e 



1 



iV m f^ oo V0 + 2Eti(l-cos ft ) + (Y} 



(Alb) 



Now if we assume that the system is infinite in d! dimensions, then we may write 
Eq. ( Alb ) in the following form 



-tTd'—d 00 

Wd ^ L ' t) = l^W £ E 



y d d 'qj dxexp j— 2 ^0 + 2^(1 — cos g;) + 



/27rt\ 2 9 N 

(— ) m . 



To obtain the last expression use has been made of the representation 



(A2) 



z Jo 



exp(—zx)dx 



and that N = L d . 

Now by rearranging it is possible to write Eq. ( |A2|) in the following form 



°° />oo 

Wd(<f>,L,t) =t E / rfxex P 



-x 40 + 2d + I — — I m' 



V A ) 



X 



— exp(2x cos q) 

L q 



d-d 1 



Here Iq(x) is the modified Bessel function. 
The use the of Poisson summation formula 



E G 



27rn 



E 



2^ 



G(g) exp(iqlL), 



(A3) 



(A4) 



(A5) 



where G(g) is a periodic function, allows us to continue the sum over the wave vector 
q = 2iin/L (n 6 [— L/2,L/2]) to the rest of the real line periodically. With the aid of 
Eq. (|A~5|) we can transform (|A4|) into 



°° />oo 

W d (0,L,t) =t E / rfxex P 



-x 4 + 2d + 



E ^il(2x) 

Z=— oo 



(A6) 



In order to investigate the low-temperature effects for the model Hamiltonian ( |2.1|) we 
use the Jacobi identity 
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^ exp ^-um 2 j = V ) E 6X ^ 



Tfl— — CO 



,7T 



-m 



m=— 00 



(A7) 



which applied to Eq. ( |A~6| ) gives 

A 



W d (<p,L,t) = XW d (<j>) + 



x 



E 



27TV2 

dx 



In 7*1/2 

m,l(d-d') u ^ 



exp 



[/ (2x)f Mac), (A8) 



where we have used the formal notations 



E Iil(2x) 

l(d-d') 



-1 d-d' 



E^(2 



x 



L I 



Z 2 — Z 2 + ■ • • + l\_ d i, 



Wd(<f>) 



2(2tt) 



j f"K fTT f _d \ 

j_ d( ii • ■ ■ / ^ u + 2 E( x - c ° s ?i) j 



'1/2 



(A9) 



the prime means that the vector with components m — l\ 



d-d' 



is omitted. 



asymptotic form for the Bessel functional 



At sufficiently low-temperature (j > 1) and large enough size [L ^> 1) we can use the 



(A10) 



7 x—v 2 /2x 



V2 



TlX 



1 9 - 32u 2 
+ 8^ + 2!(8a;) 2 + 



in order to get after substitution in Eq. (|Ala|) 

W,(d-l)/2 / K ^ 

l = x ^ + jhmr^ E 

\* n > m,l(d-d>) 



,1 2(fA) 2 m 2 + jL 2 Z 2 V 1 



2 \ V 2 

)V2 ( (A J m 2 + L 2 Z 2 



(All) 



The Watson type integral W d (4>) (see Eq. (|A"9|)) has been studied in considerable details;! 2 -^ 
for 1 < d < 3, it can be approximated by (0 <C 1), 



W d (0) « W d (0) - (4vr)-( d+1 )/ 2 



1 - d 



)/2 



(A12) 



which leads one to conclude that at zero temperature the system exhibits a phase transition 
driven by the parameter A at the quantum critical point: 



A, 



1 



W d (0)" 



(A13) 



Finally, substituting ( |A12|) in ([All] ) we obtain Eq. (f4.1|). Eq. Q3.2|) is obtained by setting 
d = d'. 
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APPENDIX B: 



In this appendix we will sketch a way to find the asymptotic behavior of the functions 
)C a (u\p,y) defined in section [TV] (see Eq. (|4.12j )). They have the following form 



fca(v\p,v)= E 



K v [2yy/r + a 2 m 2 



i(p-i) y\Jl + a 2 m 2 



y>o, 



where 



By the use of the integral representation of the modified Bessel function 

1 (z\ v l 2 

2 It, 



K v (2yfzt) = K_ v (2\[zi 



-u-l e -tx-,/* dx 



and the Jacobi identity for a p-dimensional lattice sum 

V- e -(J 2 +a 2 m 2 )t = 1 (ny r2 



m,l(p— 1) 

we may write flBlp as 



a \ t 



7 m,l(p-l) 



7r 2 (i 2 +m 2 /a 2 )/t 



(Bla) 

(Bib) 
(B2) 

(B3) 



2a 



-r 



+ 



T^-p/2 



2a 



-y 



-2u 



j-x(l 2 +rn 2 /a 2 ) 



7T\ 



p/2 



a I — 

,2 



. (B4) 



fOO , 9/0 

Jo ,, „ 

_m,i(p-l) 

Let us notice that the two terms in the square brackets in the last equality cannot be 
integrated separately, since they diverge. Nevertheless, in order to encounter this divergence, 
we can transform further ( |B4[ ) by adding and subtracting the unity from exp(— xy 2 /n 2 ), 
which enables us to write down (after some algebra) the result 



Ka{v\p,y) 



71 



p/2 



2a 



7T 



2v-p/2 i 



7T 



2i/- 



(§-!/) 



2a 



,2;/ 



E 



m,l(p— 1) 



n? 



2 \ ^-p/ 2 ' 



where 



/■oo 

C a (p|z/) = hm.y dxx 



E e 

m,l(p-l) 



-x(l 2 +m 2 /a 2 ) 



1) 

X J 



p/2 



(B5) 



(B6a) 



i im . E 



r 



8^0 



■»(P-D (* 2 + 



2 \v-p/2 



oo r 

dmd p ~H- 

-oo 



£ — ZA 



(< 2 + ^ 



iz-p/2 



(B6b) 
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is the Madelung-type constant and r[a,x] is the incomplete gamma function. 

We see from Eq. (|B5| ) that the shift of the critical quantum parameter is given by the 
Madelung type constant ( |B6| ) instead of the sum in Eq. (|4.2j). Indeed it is possible to show 
that these two representations are equivalent. This may be done, following Ref. |30|, by 
starting from the Jacobi identity Eq. (|B3|), where we multiply the two sides by §pI 2 ~ v ~ x and 
integrating over 5 to obtain the key equation 



CaipW) 



E 

m,l(p-l) 



r 



— V 



sp/ 



2-v 



I 2 + 



2 \p/2-v 



+arc 



p/2-2v 



E 



T 2 (l 2 +a 2 m 2 
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Finally form Eq. (p7|) we see easily that the integration constant C a (p\v) may be written 
in two different forms. In the first case we take the limit 5 —>■ oo and obtain 



C a (p\v)=a^/ 2 - 2 »T(v) £ 
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(B8) 



In the other case we take the limit S — > 0, and then both first and last terms in the rhs of 
Eq. (^) yields Eq. QB§). 

Using a similar procedure we find, for the functions fC a (u\p, y) defined in ( |4.13b|) , the 
following expression 
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Here the Madelung type constant is given by 
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Eqs. ( |B~5| ) and ( |B~9| ) are slight generalizations (for the anisotropic case a ^ 1) of the result 
obtained in Ref. |30] from one side, and get in touch with the Watson type sums proposed 
earlier in Ref. [13] from the other (see also Ref. [H5|) . 
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If we put in Eqs. ([B5|) or ( |B9|) d = 2, d! = and a = 1 we obtain the identity 



/f + / 2 2 y V27 V2 



+27T e1^= 1 r=J- ( Bn ) 
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FIGURES 




FIG. 1. The dependence of the universal constant yo upon the dimensionality d. The constant 
= 0.962424... is obtained for the two-dimensional system (see Eq. ( |3.12cj )) 
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FIG. 2. The effects of the finite-size geometry on the bulk behavior of cf) 1 ' 2 for the two dimen- 
sional case at A = A c . The superscript d! in Yf = Lcf) 1 / 2 and Yf = ^f^l 2 indicates the number 
of infinite dimensions in the system. The scaling variable a = The universal numbers are 
9 = 0.962424... (see Eq. (gl2^) ) and Q = 1.511955... 
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FIG. 3. The same as in FIG. ^ but for A = XtL and a = j^. The universal numbers are: 
3 = 7.061132..., St = 6.028966... and S L = 4.317795.... 
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